CHAPTER 9 FIRST-ORDER DIFFERENTIAL EQUATIONS 


9.1 SOLUTIONS, SLOPE FIELDS AND EULER'S METHOD 

1. y' = x + y => slope of 0 for the line y = —x. 

For x, y>0, y' = x + y =>■ slope > 0 in Quadrant I. 

For x, y<0, y' = x + y=> slope < 0 in Quadrant III. 

For |y| > |x|, y>0, x<0, y' = x + y=> slope > 0 in 
Quadrant II above y = —x. 

For |y| < |x|, y>0, x<0, y' = x + y => slope < 0 in 
Quadrant II below y = —x. 

For |y| < |x|, x>0, y<0, y' = x + y => slope > 0 in 
Quadrant IV above y = —x. 

For |y| > |x|, x>0, y<0, y' = x + y=í> slope < 0 in 
Quadrant IV below y = — x. 

All of the conditions are seen in slope field (d). 

2. y' = y + 1 => slope is constant for a given valué of y, slope 
is 0 for y = — 1, slope is positive for y > 1 and negative for 
y < — 1. These characteristics are evident in slope field (c). 


3. y' = — ^ =>■ slope = 1 on y = — x and —1 on y = x. 
y 1 = — * =>■ slope = 0 on the y-axis, excluding (0, 0), 
and is undefined on the x-axis. Slopes are positive for 
x > 0, y < 0 and x < 0, y > 0 (Quadrants II and IV), 
otherwise negative. Field (a) is consistent with these 
condi tions. 


4. y' = y 2 — x 2 => slope is 0 for y = x and for y = —x. 
For |y| > |x| slope is positive and for |y| < |x¡ slope is 
negative. Field (b) has these characteristics. 






Copyright © 2010 Pearson Education, Inc. Publishing as Addison-Wesley. 



538 Chapter 9 First-Order Differential Equations 


5. 


y 


iimn 

1 !!!l 

! 


m% 

i! 

t 

ÍÍÍÍT 

i 


6. 



7. y = —1 + / i X (t-y(t))dt^g=x-y(x);y(l) = -l + £ (t - y(t))dt = -1; g = x - y, y(l) = -1 

8 - y = /> => I = i; y(l) = {dt = 0 ; g = {, y(l) = 0 

9. y = 2 - J o X (l +y(t))sintdt =» g = -(1 + y(x))sinx; y(0) = 2 - £(1 + y(t))sintdt = 2; g = -(1 +y)sinx, 

y (0) = 2 

10. y = 1 + £ y(t) dt => | = y(x); y(0) = 1 + £ y(t) dt = 1; g = y, y(0) = 1 

11. Yi = Yo + (l - g) dx = -1 + (l - (.5) = —0.25, 

ya = yi + (i - g) dx = -0.25 + (1 - =§¥) (.5) = 0.3, 

y 3 = y 2 + (l - g) dx = 0.3 + (1 - f ) (.5) = 0.75; 

g + (i) y = 1 P(x) = i , Q(x) = 1 => f P(x) dx = /{ dx = ln |x| = ln x, x > 0 v(x) = e lnx = x 
=> y = 1 fx ■ 1 dx = i (f + c) ; x = 2, y = -1 => -1 = 1 + § C = -4 => y = § - \ 

=> y(3.5) = = - y 5 « 0.6071 

12. yi = y 0 + x 0 (1 - y 0 ) dx = 0 + 1(1 - 0)(.2) = .2, 

y 2 = yi + X! (1 - yi) dx = .2 + 1.2(1 - ,2)(.2) = .392, 

y 3 = y 2 + X 2 (1 - y 2 ) dx = .392 + 1.4(1 - ,392)(.2) = .5622; 

= x dx => —ln|l — y| = y+C;x=l, y = 0 => -ln 1 = 1 + C => C = - \ => ln|l-y| = -^ + i 
=> y = 1 - e (1 ~ x2)/2 y(1.6) « .5416 

13- yi = y 0 + (2x 0 y 0 + 2y 0 ) dx = 3 + [2(0)(3) + 2(3)](.2) = 4.2, 

y 2 = Yl + (2x iyi + 2 yi ) dx = 4.2 + [2(.2)(4.2) + 2(4.2)](.2) = 6.216, 

y 3 = y 2 + (2x 2 y 2 + 2y 2 ) dx = 6.216 + [2(.4)(6.216) + 2(6.216)](.2) = 9.6969; 

g = 2y(x +1) => f = 2(x + 1) dx =4- ln |y| = (x + l) 2 + C; x = 0, y = 3 =4> ln 3 = 1 + C =4> C = ln 3 - 1 

=> ln y = (x + l) 2 + ln 3 - 1 =4> y = e ( x+ 0 2 + 1 " 3 - 1 = e ln3 e x2+2x = 3e x ( x+2) => y(.6) « 14.2765 

14. y x = y 0 + y§(l + 2x 0 ) dx = 1 + 1 2 [1 + 2(-l)](.5) = .5, 

y 2 = yi + y?(l + 2x0 dx = .5 + (,5) 2 [1 + 2(-.5)](.5) = .5, 

y 3 = y 2 + y|(l + 2x 2 ) dx = .5 + (,5) 2 [1 + 2(0)](.5) = .625; 

g = (1 +2x)dx =4> - I = x + x 2 +C;x = -l,y = 1 =4> -1 = -1 +(-l) 2 + C => C = -l => 1 = 1 - x- 

^ y = i - x - x 2 ^ y(- 5 ) = i—. 5 —(.5) 2 = 4 
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15. y x = y 0 + 2x 0 e x ° dx = 2 + 2(0)(.l) = 2, 

y 2 = Yl + 2xie x5 dx = 2 + 2(.l) e l2 (.l) = 2.0202, 
y 3 = y 2 + 2 x 2 e x 2 dx = 2.0202 + 2(.2) e- 22 (. 1) = 2.0618, 

dy = 2xe x2 dx => y = e x ' + C; y(0) = 2=+2=l+C=í>C=l => y = e x2 + 1 => y(.3) = e 3 ' + 1 « 2.0942 

16. yi = yo + (yo e x °) dx = 2 + (2 • e°) (.5) = 3, 

y 2 = yi + (yi e Xl ) dx = 3 + (3 • e 05 ) (.5) = 5.47308, 

y 3 = y 2 + (y 2 e X2 ) dx = 5.47308 + (5.47308 • e 10 ) (.5) = 12.9118, 

= y e x => ^ = e x dx =$■ ln|y| = e x + C; x = 0, y = 2 =$■ ln 2 = 1 + C =>• C = ln 2 — 1 =+ ln|y| = e x + ln 2 — 1 

=+ y = 2e eX - 1 =>• y( 1.5) = 2e e ' 5 - 1 « 65.0292 


17. yi = 1 + 1 (.2) = 1.2, 

y 2 = 1.2 + (1.2)(.2)= 1.44, 
y 3 = 1.44 +(1.44)(.2) = 1.728, 
y 4 = 1.728 + (1.728)(.2) = 2.0736, 
y 5 = 2.0736 + (2.0736)(.2) = 2.48832; 

^ = dx => lny = x + Ci => y = Ce x ; y(0) =1 =+ 1 = Ce 0 =+ C = 1 => y = e x => y(l) = e « 2.7183 

18. yi =2 + (f) (.2) = 2.4, 

y 2 = 2.4 + (y=|) (.2) = 2.8, 
y 3 = 2.8 + (f|) (.2) = 3.2, 
y 4 = 3.2+ (f|) (.2) = 3.6, 
y 5 = 3.6 + (jJ) (.2) = 4; 

f = f => ln y = ln x + C => y = kx; y(l) = 2=+2 = k=>y = 2x=>- y(2) = 4 


19. yi = —1 + 
y 2 = --5 + 


(-D 2 


L vd J 

(—. 5 )' 


vdíj 

y 3 = -.39794 + 
y 4 = -.34195 + 


05) = -.5, 

(.5) = -.39794, 

(.5) = -.34195, 
(.5) = -.30497, 


(—. 39794) 2 
y/2 

(—. 34195) 2 


y/25 

y 5 = -.27812, y 6 = -.25745, y 7 = -.24088, y 8 = -.2272; 

| = ^ =► -i=2 v ^ + C;y(l) = -l => 1=2 + C => C = 


1 ^ y_ í-lyfi ^ y(5)~ j_ 2a/ 5 


.2880 


20. yi = 1 + (0 • sin 1) (i) = 1, 

y 2 = 1 + (I . sin l) (i) = 1.09350, 
y 3 = 1.09350 + (| • sin 1.09350) (§) = 1.29089, 

y 4 = 1.29089 + (| • sin 1.29089) Q) = 1.61125, 

y 5 = 1.61125 + (| • sin 1.61125) Q) = 2.05533, 

y G = 2.05533 + (f • sin2.05533) (§) = 2.54694; 

y' = xsiny => cscydy = xdx => —ln|cscy + coty| = ^x 2 + C =+ cscy + coty = e _ 5 x2+c = Ce^5 x2 
=> 1 ^5° sy = Ce~í x " => cot( y ) = Ce _ J x2 ; y(0) = 1 => cot()) = Ce 0 = C =+ cot ( 2 ) = cot(i)e _ i x2 

=> y = 2cot _l (cot(j)e _ í x '), y(2) = 2cot _I (cot(|)e~ 2 ) = 2.65591 
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21. y = -1 - X+ (1 +x 0 + yo)e x x ° => y(x 0 ) = -1 - x 0 + (1 + x 0 + y 0 )e x ° x « =-1 - x 0 + (1 + x 0 + y 0 )(l) = yo 
% = -1 + (1 +x 0 + y 0 )e x “ x ° =>• y = ~1 -x + (1 +x 0 + y 0 )e x " Xo = ^- x=»^=x + y 

22. y' = f(x), y(x 0 ) = y 0 =>■ y = f* f(t)dt + C, y(x 0 ) = f(t)dt + C = C=^C = y 0 =>y = fjm +y 0 


23-34. Example CAS commands: 

Maple : 

ode := diff( y(x), x ) = y(x); 
icA := LO, 1]; 
icB := [0. 2]; 
icC := [0,-1]; 

DEplot( ode, y(x), x=0..2, [icA,icB,icC], arrows=slim, linecolor=blue, title="#23 (Section 9.1)"); 
Mathematica : 

To plot vector fields, you must begin by loading a graphics package. 

«Graphics'PlotField' 

To control lengths and appearance of vectors, select the Help browser, type PlotVectorField and select Go. 
Clear[x, y, f] 
yprime = y (2 — y); 

pv = PlotVectorField[{ 1, yprime}, ¡x, —5, 5}, {y, —4, 6], Axes —> Trae, AxesLabel —> {x, y}]; 

To draw solution curves with Mathematica, you must first solve the differential equation. This will be done with 
the DSolve command. The y[x] and x at the end of the command specify the dependent and independent variables. 
The command will not work unless the y in the differential equation is referenced as y[x]. 
equation = y'[x] == y[x] (2 — y[x]) ; 
initcond = y[a] == b; 

sois = DSolve[{equation, initcond], y[x], x] 

vals = {{0, 1/2], [0,3/2], [0,2], {0,3}} 

f[{a_, b_¡] = sols[[ 1, 1,2]]; 

solnset = Map[f, vals] 

ps = Plot [Evalúate [solnset, [x, —5, 5}]; 

Show[pv, ps, PlotRange —> {—4,6}]; 

The code for problems such as 31 & 32 is similar for the direction field, but the analytical Solutions involve 
complicated inverse functions, so the numerical solver NDSolve is used. Note that a domain interval is 
specified. 

equation = y'[x] == Cos[2x — y[x]] ; 
initcond = y [0] == 2; 

sol = NDSolve [{equation, initcond], y[x], {x, 0, 5}] 
ps = Plot[Evaluate[y[x]/.sol, {x, 0, 5}]; 

N[y[x] /. sol/.x -► 2] 

Show[pv, ps, PlotRange —> {0,5}]; 

Solutions for 34 can be found one at a time and plots named and shown together. No direction fields here. 

For 34, the direction field code is similar, but the solution is found implicitly using integrations. The plot 
requires loading another special graphics package. 

«Graphics' ImplicitPlot' 

Clear[x,y] 

solution[c_] = lntegrate[2 (y — 1), y] == Integrate[3x 2 + 4x + 2, x] + c 

valúes = {—6, —4, —2, 0, 2, 4, 6}; 

solns = Maplsolution, valúes]; 

ps = ImplicitPlot[solns, {x, —3, 3], {y, —3, 3}] 

Show[pv, ps] 


Copyright © 2010 Pearson Education, Inc. Publishing as Addison-Wesley. 



Section 9.1 Solutins, Slope Fields and Euler's Method 541 



27. 28. 



35. ^ = 2xe x \y(0) = 2 => y n+ i = y n + 2x n e x »dx = y n + 2x n e x »(0.1) = y n + 0.2x n e x » 

On a TI-84 calculator home screen, type the following commands: 

2 STO > y: 0 STO > x: y (enter) 

y + 0.2*x*e A (x A 2) STO > y: x + 0.1 STO > x: y (enter, 10 times) 

The last valué displayed gives yEuier(l) ~ 3.45835 

The exact solution: dy = 2xe x2 dx =£• y — e x2 I - C; y(0) — 2 = e° + C => C — 1 => y = 1 + e x2 
yexact (1) = l + e« 3.71828 

36. 37 = 2y 2 (x - 1), y(2) = =>• y n +i = y n + 2y 2 (x n - l)dx = y n + 0.2y 2 (x n - 1) 

On a Tl-84 calculator home screen, type the following commands: 

—0.5 STO > y: 2 STO > x: y (enter) 

y + 0.2*y 2 (x — 1) STO > y: x + 0.1 STO > x: y (enter, 10 times) 

The last valué displayed gives yEuier(2) ~ —0.19285 

The exact solution: ^ = 2y 2 (x — 1) => p) = (2x — 2)dx =>—l=x 2 — 2 x + C=> 3 = _ x 2 4 - 2x + C 
y(2) = -3 => rr /2 = -(2 ) 2 + 2(2) + C = C => C = —2 3 = —x 2 + 2x — 2 => y = _ x2+ ‘ 2x _ 2 

“ — ( 3) 2 + 2 ( 3 ) — 2 = _0 ' 2 
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37. s = y > 0, y(0) = l =* y„+i = y„ + ^dx = y n + ^(0.1) = y n + 0.1^ 

On a TI-84 calculator honre screen, type the following commands: 

1 STO > y: 0 STO > x: y (enter) 

y + 0.1*(^/x /y) STO > y: x + 0.1 STO > x: y (enter, 10 times) 

The last valué displayed gives yEuier(l) ~ 1.5000 

The exact solution: dy = ^dx y dy = ^/x dx => \ = |x 3/2 + C; = % — \ = §(0) 2/2 + C => C = \ 

=► i = |x 3/2 + ¡ => y = -/fx^TT =► y exact (1) = ^(l) 3/2 + l « 1.5275 

38. | = 1+ y 2 , y(0) = 0 => y n+1 = y n + (1 + y 2 )dx = y„ + (1 + y 2 )(0.1) = y„ + 0.1(1 + y 2 ) 

On a Tl-84 calculator home screen, type the following commands: 

0 STO > y: 0 STO > x: y (enter) 

y + 0.1 *( 1 + y 2 ) STO > y: x + 0.1 STO > x: y (enter, 10 times) 

The last valué displayed gives yEuier(l) ~ 1.3964 

The exact solution: dy = (1 + y 2 )dx => | 2 = dx => tan _1 y = x + C; tan _1 y(0) = tan^'0 = 0 = 0 + C=>C = 0 

=> tan~'y = x => y = tan x => y ex act(l) = tan 1 « 1.5574 

39. Example CAS commands: 

Maple : 

ode := diff( y(x), x ) = x + y(x);ic := y(0)=-7/10; 
xO := -4;xl := 4;y0 := -4; yl := 4; 
b := 1; 

P1 := DEplot( ode, y(x), x=x0..xl, y=y0..yl, arrows=thin, title="#39(a) (Section 9.1)"): 

Pl; 

Ygen := unapply( rhs(dsolve( ode, y(x))), x,_Cl ); # (b) 

P2 := seq( plot( Ygen(x.c), x=x0..xl, y=y0..yl, color=blue ), c=-2..2 ): # (c) 

display( [Pl,P2], title="#39(c) (Section 9.1)"); 

CC := solve( Ygen(0,C)=rhs(ic), C ); # (d) 

Ypart := Ygen(x.CC); 

P3 := plot( Ypart, x=0..b, title="#39(d) (Section 9.1)"): 

P3; 

euler4 := dsolve( {ode,ic}, numeric, method=classical[foreuler], stepsize=(xl-x0)/4 ): # (e) 

P4 := odeplot( euler4, [x,y(x)], x=0..b, numpoints=4, color=blue ): 
display( [P3,P4], title="#39(e) (Section 9.1)"); 

euler8 := dsolve( {ode,ic}, numeric, method=classical[foreuler], stepsize=(xl-x0)/8 ): # (f) 

P5 := odeplot( euler8, [x,y(x)], x=0..b, numpoints=8, color=green ): 

eulei'16 := dsolve( {ode.ic}, numeric, method=classical[foreuler], stepsize=(xl-x0)/16 ): 

P6 := odeplot( eulerló, [x,y(x)j, x=0..b, numpoints=16, color=pink ): 

euler32 := dsolve( {ode.ic}, numeric, method=classical[foreuler], stepsize=(xl-xO)/32 ): 

P7 := odeplot( euler32, [x,y(x)j, x=0..b, numpoints=32, color=cyan ): 
display( [P3,P4,P5,P6,P7], title="#39(f) (Section 9.1)"); 

< < N | h | 'percent error' >, # (g) 

< 4 | (xl-x0)/ 4 | evalf[5]( abs(l-eval(y(x),euler4(b))/eval(Ypart,x=b))*100 ) >, 

< 8 | (xl-x0)/ 8 | evalf[5]( abs(l-eval(y(x),euler8(b))/eval(Ypart,x=b))*100 ) >, 

< 16 | (xl-x0)/16 | evalf[5]( abs(l-eval(y(x),eulerl6(b))/eval(Ypart,x=b))*100 ) >, 

< 32 | (xl-x0)/32 | evalf[5]( abs(l-eval(y(x),euler32(b))/eval(Ypart,x=b))*100 ) > >; 
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39-42. Example CAS commands: 

Mathematica : (assigned functions, step sizes, and valúes for initial conditions may vary) 

Problems 39 - 42 involve use of code from Problems 23 - 34 together with the above code for Euler's method. 


9.2 FIRST-ORDER LINEAR DIFFERENTIAL EQUATIONS 

1. xg+y = e x ^g + (i)y=f,P(x)=I,Q(x)=f 

/ P(x) dx = fí dx = ln |x| = ln x, x > 0 => v(x) = e/ p(x)dx = e lnx = x 
y=é)/ v WQ« dx= x/x(f)dx=i(e x + C) = ^,x>0 


2. e x g + 2e x y = 1 =► g + 2y = e~ x , P(x) = 2, Q(x) = e~ x 
/P(x) dx = J 2 dx = 2x =$■ v(x) = e/ p(x) dx = e 2x 

y = ¿r Je 2x • e- x dx = ^ /e x dx = ^ (e x + C) = e~ x + Ce- 2x 

3. xy' + 3y = ^ , x > 0 =* g + (5) y = , P(x) = 3 , Q(x) = si V 

/ | dx = 3 ln |x| = ln x 3 , x > 0 => v(x) = e lnxi = x 3 

y = F jx 3 (^) dx = i /sin x dx = ¿ (- eos x + C) = . x > 0 

4. y' + (tan x) y = eos 2 x, — |<x<f=»^-l- (tan x) y = eos 2 x, P(x) = tan x, Q(x) = eos 2 x 

/tan x dx = / dx = — ln |cos x| = ln (eos x)~ x , — f < x < | =>■ v(x) = e lnícosx) = (eos x)” 1 

y = (cos 1 xr i /(eos x) _1 - eos 2 x dx = (eos x)/eos x dx = (eos x)(sin x + C) = sin x eos x + C eos x 


dy 

dx 

'2 


2y=l — ^,x>0=> 


+ (i) y = x - P(x > = x > QW = x 


/= dx = 2 ln |x| = ln x 2 , x > 0 =>• v(x) = e lnx2 = x 2 


y = 


r /x 2 (I-¿)dx=¿/(x-l)dx=i(f-x + C ) = I-I + § ) x>0 


6 . 


(i + x)y' + y = v^^l + ( T ri)y 


jA 

1+ x 


,P(X) 


1 

1+x 


, Q(X) = 


jAL 

1+x 


/dx = ln (1 + x), since x > 0 =>■ v(x) = e ln(1+x) = 1 


y = 


1+x 


/(l + X) (4) dx = tC / A * = (rC) (i x"' 2 + C) 


2x 3 / 2 , C 

3(1+x) ' t ' 1+x 


7. S - \ y = l eX/2 => P(x) = - \ , Q(x) = ¿ e x/2 =+ /P(x) dx 
=> y = ^ /e-x /2 (i e */2) dx = e x/2 / \ dx = e x/2 (| x + C) 


-|x v(x) = e x/2 
| xe x/2 + Ce x/2 


8. g + 2y = 2xe- 2x => P(x) = 2, Q(x) = 2xe- 2x => /P(x) dx = /2 dx = 2x => v(x) = e 2x 

=>■ y = i/e 2x (2xe~ 2x ) dx = ^/ 7x dx = e~ 2x (x 2 + C) = x 2 e~ 2x + Ce 2x 

9. g - (i) y = 2 ln x =>• P(x) = - i , Q(x) = 2 ln x => / P(x) dx = -/i dx = - ln x, x > 0 

=+ v(x) = e~ lnx = i =>■ y = x/ (1) (2 ln x) dx = x [(ln x) 2 + C] = x (ln x) 2 + Cx 
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10. g + (?) y = , X > 0 ^ P(x) = \ , Q(x) = ■ => fp(x) dx = /? dx = 2 ln |x| = ln x 2 , x > 0 

=> v (x) = e lnx2 = x 2 => y = 4 Jx 2 (^) dx = i/eos x dx = j¡ (sin x + C) = 

H. I + (f4r) S = ^3 =► P(t) = , Q(t) - (T^e =* /P(t) dt = / i 4 T dt = 41n|t—1| = ln(t - l) 4 

=> v(t) = e ln ^ 1)l = (t - l) 4 => s = ^4 f(t- l) 4 [^rryr] dt = /(t 2 - 1) dt 

1 {$ t C 

~ (t-l) 4 ^ 3 ~ 1 d” CJ = 3(t- ly 1 — (t- l) 4 + (t-1 f 

12. (1+1)1+28 = 3(1+1)+^ =s- | + (f/r) 8 = 3 + ^ =► P(t)= r ^ T ,Q(t) = 3 + (t+l)- 3 

=> /P(t) dt = dt = 2 ln |t + 1| = ln(t + l) 2 => v(t) = e ln(t+1 > 2 = (t + l) 2 

s = /(t + l) 2 [3 + (t + 1) 3 ] dt = (f^yi /[3(t + l) 2 + (t + 1) 4 ] dt 

= (T+T)i [(t + l) 3 + ln |t + 11 + C] = (t + 1) + (t + 1) 2 ln(t + 1) + )2 , t > — 1 

13. ^ + (cot 0)r = sec 0 =>■ P (6) = cot 9, Q (9) = sec 9 => /P(0) d9 = /cot 9 d9 = ln |sin 9\ => v(0) = e 1 ” 18 '”" 1 

= sin 9 because 0 < 0 < | =>■ r = —L / ( s i n g)( sec $) d9 — —L / t an Q d9 = J-g (ln |sec 9\ + C) 

= (esc 9) (ln |sec 9\ + C) 

14. tan 9 % + r = sin 2 0 => af + ^ = =4> ^ + (cot 9) r = sin 9 eos 9 => P (9) = cot 9, Q(0) = sin 9 eos 9 

=> /P(0) d0 = /cot 0 d0 = ln |sin 9\ = ln(sin 9) since 0 < 9 < | => v(0) = e ln(smfl ) = sin 9 

=► r = ¡ib/( sin 0 )( sin 0 cos d0 = üb / sin2 6 cos 0 d0 = (¡ib) (bn + c) = H 1 + £- e 

15. | + 2y = 3 =S> P(t) = 2, Q(t) = 3 => f P(t) dt = / 2 dt = 2t => v(t) = e 2t => y = ¿ /3e 2t dt 
= b(| e2t + C);y(0)=l ^|+C=1^C = -1^ y = f - | e" 2 * 

16. | + f = t 2 =>• P(t) = ? , Q(t) = t 2 => /P(t) dt = 2 ln |t| => v(t) = e ln ' 2 = t 2 => y = £/(t 2 ) (t 2 ) dt 

= ? f t4 dt = ? (i + c ) = V + § ; y( 2 )= 1 f + í = 1 ^ c = _ t y = V — 5^ 

17. I + 0) y = ^ =► p ( 0 ) = s > W) = ^ =+ /P(0) d0 = ln |0| => v(0) = e 1 "' 8 ' = |0| 

=> y=//|0| (^) d0= l/e (sai) d0 for0^O => y = i/sin 0 d0 = i (- cos 9 + C) 

= - i cos 0 + f ; y (|) = i =+ C = | ^y = -l C os0+^ 

18. — (|) y = 0 2 sec 9 tan 9 => P(0) = — |, Q(0) = 0 2 sec 9 tan 0 => /P(0) d0 = —2 ln |0| => v(0) = e- 21 ” 1 " 1 
= 0 -2 =>■ y = g4j/(0 -2 ) (0 2 sec 0 tan 0) d0 = 0 2 /sec 0 tan 0 d0 = 0 2 (sec 0 + C) = 0 2 sec 0 + C0 2 ; 

y (f) = 2 =► 2 = (f) (2) + C (/) =► C = - 2 =* y = 0 2 sec 0 + (« - 2) 0 2 


19. (x+l)^-2(x 2 +x)y= x ^ T =* I-2[^ijy= ( ^ T7 ^ | - 2xy = ^ => P(x) =-2x 
Q(X) = ^ =+ /P(x) dx = /—2x dx = -x 2 ^ v(x) = e-* 2 => y = ^/e^ 2 [^] dx 
= eX 7üW dx = ex2 Í^F+C] =-^ I +Ce x2 ;y(0) = 5 ^ ^^+0 = 5 ^ -1+C = 5 


P(x) = —2x, 


Q(x) = 


=► C = 6 ^ y = 6e x " — 4 zt 
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20. g + xy = x =X P(x) = x, Q(x) = x => f P(x) dx = fx dx = f => v(x) = e x2 / 2 =>• y = \ f e x¡/2 - x dx 

= ^ (e xV2 + C) = 1 + ;y(0) = -6 ^ 1 +C = -6 ^ C = -7 ^ y = 1 - ^ 

21. ^ - ky = 0 => P(t) = -k, Q(t) = 0 => f P(t) dt = f k dt = -kt => v(t) = e~ kt 

=> y = ¿hs f (e” kt ) (0) dt = e kt (0 + C) = Ce kt ; y(0) = y 0 =xC = y 0 =xy = y 0 e kt 

22. (a) | + k u = 0 ^ P(t) = £, Q(t) = 0 => f P(t) dt = dt = £ t = £ =► u(t) = e k '/ m 

=> y = etk fe kt / m - o dt = ggjs; u(0) = u 0 =x ¿¿§ 7 = = u 0 =x C = u 0 =>• u = u 0 e~k/m)t 
(b) ^ = — — u — = — — dt => l nu = ——t+C U = e - ( k / m k +c => u = e~k/ m ) 1 . e c . Lete c = Ci. 

v/ dtmum m 1 

Then u = g^ • Ci and u(0) = u 0 = g^ • Ci = Ci. So u = u 0 e"^" 1 ) 1 

23. x J'^ dx = x (ln |x| + C) = x ln |x| + Cx => (b) is correct 


24. 


1 

eos X 


I eos x dx = —— (sin x + C) = tan x H—— 

J eos x v J eos x 


=> (b) is correct 


25. Steady State = k and we want i = | (k) =>■ | (jQ = X (1 _ e Rt / L ) I — \ _ e Rt /L — i = —e Rt / L 

=> ln | = — => — ^ ln i = t => t = g ln 2 sec 

26. (a) | £ i = 0 =X i di = - l dt =>• ln i = - f + Cj =X i 

=> i = Ie~ Rt/L amp 

(b) \\ = I e - Rt / L =>• e - Rt / L = 1 =>• - f = ln i = - ln 2 =X t 

(C) t = g =>. i = I e (- R t/L)(L/R) _ J g-t am p 

27. (a) t — 3 j k =>- i = g (1 — é R/L )( 3L / R j) = X (l — e~ 3 ) « 0.9502 k amp, or about 95% of the steady State valué 

(b) t = ^ => i = k (l — 6 (-m.)(2l/r)^ _ v _ e - 2 ^ ^ 0.8647 k amp, or about 86% of the steady State valué 


= e Cl e~ Rt/L = Ce- Rt / L ; i(0) = I 1 = C 
= ^ ln 2 sec 


28. (a) 


(c) 


di 

dt 


V 

L 


P(t) = f , Q(t) = 


/ P(t) d. = / 


R dt = — 

L UL L 


v(t) = e 


,Rt/L 


1 = 


1 

e Rt/L 


j> ,L (B dt = * [| © + C] = | + Ce->'» 


(b) i(0) = 0 ^ £+C = 0=>C = -£ => i = 


i _ V _ V p—Rt/L 


R 


R 


i=í 


gj = 0 => || + R i = 0 + (c) (k) = E i — ^ is a solution of Eq. (11); i = Ce (R/L)t 


29. y' — y = —y 2 ; we have n = 2, so let u = y 1-2 = y -1 . Then y = u -1 and ^ = — ly _2 gj => gf = — y 2 |k 
=> _ u -2 ^ — u _1 = —u -2 =>• ^ + u = 1. With el dx = e x as the integrating factor, we have 
eX (s + u ) = ®( eXu ) = eX - Integrating, we get e x u = e x + C =>• u = 1 + £ = I => y = 


30. 


y' - y = xy 2 ; we have n = 2, so let u = y” 1 . Then y = u~' and g = -y" 2 ^ => g = -y 2 g = -u" 2 ^. 
Substituting:—u -2 — u _1 = xu~ 2 => + u = —x. Using e/ dx = e x as an integrating factor: 


+u ) = ¿(e x u) = -xe" =$■ e x u = e x (l - x) + C 


u = 




- xe x + C 


31. xy' + y = y~ 
= 3y 2 


=> y' + (I)y = Q)y~ 2 . Let u = y 1_ (“2) _ y¿ y = u uj an( j _ u - 

y' = ‘¿ = (IKsiXr 2 ) = (jXsH^Thuswehave 


du 

dx 


1? 2dy 

dx 


2/3 
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( 5 ) (gj) (u 2/3 ) + (i)u 1/3 = (j)u 2/3 => gj + (j)u = (|) 1. The integrating factor, v(x), is 
eix dx = e 31nx = e lnx3 = x 3 . Thus j ¡ ;(x 3 u) = (f)x 3 = 3x 2 => x 3 u = x 3 +C=t>u=l + L(=y 3 

=>y=(i + §) 1/3 

32. x 2 y' + 2xy = y 3 ^ y'+ (|)y = (¿)y 3 . P(x) = (=),Q(x) = (¿), n = 3. Let u = y '- 3 = y" 2 . 

Substituting gives ^ + (—2)( = )u = —2(¿) => gj + (g^)u = g(r- Let the integrating factor, v(x), be 
e/N) dx = e lnx 4 = x~ 4 . Thus ¿(x -4 u) = —2x -6 => x " 4 u = |x -5 + C => u = L;+Cx 4 = y -2 
^y=(it+cx 4 )- 1/2 

9.3 APPLICATIONS 

1. Note that the total mass is 66 + 7 = 73 kg, therefore, v = voe~^ m ^ => v = 9 e ~ 3 - 9t / 73 

(a) s(t) = /9e -3 - 9t / 73 dt = -2i|Q e -3.9t/73 + c 

Since s(0) = 0 we have C = ^f 9 and lims(t) = lim 2 ¡ 90 (l - e~ 3 ' 9t/73 ) = « 168.5 

The cyclist will coast about 168.5 meters. 

(b) 1 = 9e ~ 3 - 9t / 73 =s> L|t = in 9 =>. t = ^2 ~ 41.13 sec 
It will take about 41.13 seconds. 

2 v = v 0 e“( k / m)t => v = 9 e-( 59 - 000 / 51 > 000 ' (XX) )t v = 9 e ~ 59 t / 5 i.ooo 

(a) s(t) = /9e —59t / 51 000 dt = -«Mx» e -59t/5i,ooo + c 

Since s(0) = 0 we have C = ^M ® 99 and lims(t) = lim l 39 ^ 99 (l - e - 59 t/ 5 i,ooo\ = «mmoo _ 7780 m 

Hoo Neo jy \ / Dy 

The ship will coast about 7780 m, or 7.78 km. 

(b) 1 = 9 e ” 59t/51 ’ 000 => NjN = | n 9 ^ t = 51 -°°° ln9 w 1899.3 sec 
It will take about 31.65 minutes. 

3. The total distance traveled = v “ ,ri => l 2 - 75 ^ 39 - 92 ) = 4 91 => k = 22.36. Therefore, the distance traveled is given by the 
function s(t) = 4.91 (l — e ^ 22 - 3 6/39.92)t) The graph shows s(t) and the data points. 



4 . ^ = coasting distance =t> (o.80)(49.90) _ 1 32 =$. k = ^ 

We know that ^ = 1.32 and L = = §§. 

Using Equation 3, we have: s(t) = ^ (l — e N /m4 ) = 1.32(1 — e _29t/33 ) « 1.32(1 — e _0 - 696t ) 
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5. y = mx => | = m => xy x2 y = 0 => y' = So for 
orthogonals: ^ ^ => y dy = — x dx => ^- + y = C 

=> x 2 + y 2 = Ci 


6 . y = ex 2 => = c =>■ x y x 4 2xy = 0 =>• x 2 y' = 2xy 

=> y' = ^• So for the orthogonals: ^ = — A 

=> 2ydy = -xdx =^y 2 = -^+C=>y = ± yj f + C, 

C > 0 


7. kx 2 + y 2 = 1 => 1 - y 2 = kx 2 => = k 

=> x 2 ( 2 y)y ';i 1 - y2)2x = 0 =s- — 2 yx 2 y' = (1 - y 2 )( 2 x) 
=> y' = ^ -o X yj7 X) = ^_ X y - ■ So for the orthogonals: 

% = => = xdx => lny - £ = ¿ + c 


8 . 2x 2 + y 2 = c 2 => 4x + 2yy' = 0 => y' = = -y. For 

orthogonals: ^ ^ =>■ ^ ^ =>■ ln y = ^ln x + C 

=> ln y = ln x 1/2 + InCi =>■ y = Ci |x | 1/2 


n —x v e V -y(e x )(— l) f, 

9. y = ce x => = c =$■ — y /V, -- = 0 

2 e (e - 1 ) 2 


=> e x y' — —ye x =» y' = —y. So for the orthogonals: 
| ± =»ydy = dx^ ¿ =x + C 

=> y 2 = 2 x + Q => y = d= yj 2 x + Ci 




y 



y 
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, i x E)y'-iny 

10 . y = e kx =$■ ln y = kx =>■ = k => ^ x 2 - = 0 

=> y' — ln y = 0 => y' = $o f or the orthogonals: 

S = ^=^y ln y d y = - xdx 

=^b 2ln y-?(y 2 ) = (-5 x2 )+ c 

=> y 2 ln y - ¿ = -x 2 + Q 



11. 2x 2 + 3y 2 = 5 and y 2 = x 3 intersect at (1, 1). Also, 2x 2 + 3y 2 = 5 => 4x + 6 yy' = 0 => y' = — g 5 =»■ y'(l, 1) = 
y 2 = x 3 => 2yiy [ = 3x 2 =>• y{ = => y{(l, 1) = f. Since y' • y[ = (—|) (|) = — 1, the curves are orthogonal. 


2 2 

12. (a) x dx + y dy = 0 => y + \ = C is the general equation 
of the family with slope y' = —For the orthogonals: 
y' = | f =*► ln y = ln x + C or y = Qx 

(where Ci = e c ) is the general equation of the 
orthogonals. 


(b) x dy - 2 y dx = 0 2 y dx = x dy => ^ = f 

=> 1 f =>. Iln y = ln x + C =>■ y = Cix 2 is 

the equation for the solutionfainily. 

^ln y — lnx = C=>i — — - =0=>y'= — 

2 J 2 y x J x 

=>• slope of orthogonals is ^ = — j- 
=> 2y dy = —x dx y 2 — - y I C ¡s the general 
equation of the orthogonals. 



13. Let y(t) = the amount of salt in the container and V(t) = the total volunte of liquid in the tank at time t. 
Then, the departure rate is ^4 (the outflow rate). 

(a) Rate entering = = 10 lb/min 

(b) Volume = V(t) = 100 gal + (5t gal - 4t gal) = (100 + t) gal 

(c) The volume at time t is (100 + t) gal. The amount of salt in the tank at time t is y lbs. So the 
concentration at any time t is 10 ^ lbs/gal. Then, the rate leaving = 1Q ^ (lbs/gal) • 4 (gal/min) 

= 7boA lbs/min 

(d) | = 10 ^h =* I + (iW+t)y = 10 ^ p (t) = ioü+í > Q(t) = 1° =* /P(t) dt = /to^r dt 

= 4 ln (100 + t) =► v(t) = e 41n(100+t > = (100 + t ) 4 =► y = < 100 ^ /(100 + t) 4 (10 dt) 

= (TMW (^^+c) =2(100+ 1 )+^=^; y(0) = 50 =* 2(100+ 0 +^ 0 ? =50 
=► c = — (150)(100) 4 =► y = 2(100 + t)-«^! =* y = 2(100 +1) - -^0 

J V 1 + 100/ 
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(e) y(25) = 2(100 + 25) - « 188.56 lbs =* concentration = ^ ~ » L5 lb/gal 

14. (a) ^ = (5 - 3) = 2 => V = 100 + 2t 

The tank is full when V = 200 = 100 + 2t => t — 50 min 

(b) Let y(t) be the amount of concéntrate in the tank at time t. 

dy _ (l Jb (r gaA _ ( y Ib \ Gga('\^,dy = 5_ 3 ( y \ _v dy , 3 _ 5 

dt y 2 gal J y min J yl00 + 2tgal J y min J dt 2 2 \50 + t/ dt 2(t + 50)' 2 

Q(t) = P(t) = | (joVt) => /P(t) dt = dt = |ln (t + 50) since t + 50 > 0 

v(t) = e J' P(t,dt = eJ ln(t + 5 °) = (t + 50) 3/2 

y(t) = + 50) 3/2 dt = (t + 50)- 3/2 [ (t + 50) 5/2 + c] =* y(t) = t + 50 + Jg^ 

Apply the initial condition (i.e., distilled water in the tank at t = 0): 

y(0) = 0 = 50 + =>■ C = —50 s/2 = 2 - y(t) = t + 50 — - m . When the tank is full at t = 50, 

y(50) = 100 — « 83.22 pounds of concéntrate. 


15. Let y be the amount of fertilizer in the tank at time t. Then rate entering = 1 ^ 


1 gal = 1 


and the 


volume in the tank at time t is V(t) =100 (gal) + [1 (gal/min) — 3 (gal/min)]t min = (100 — 2t) gal. Henee 


rate out = I 

=► P(t) = 


y ) 

100 -2tJ 


3 = T0Ü=2t lbs/mÍn 


dy 


= 0 - 


3y 


100-21 > Q® 1 

= (100 - 2 t ) -3 / 2 


dt — V a 100 — 2t 

Jp(t) dt = / > 


lbs/min 


dy 


_ 3 ln(100—2t) 
100-21 Ul — -2 


' dt + ( loo —2t) y — 1 

v (t) = e ( 31n(100 2t))2 


y doo-2tr 3 / 2 
= (100 - 2t) + C(100 - 2t) 3 / 2 ; y(0) = 0 


f (100 - 2 t )~ 3 / 2 dt = (100 - 2 t)“ 


3/2 


— 2(100 — 2 tj ~ 1 / 2 


=> C = -(100r 1/2 = - i => y = (100 - 2t) - uuu 7n ¿u ~" . Let 1 = 0 => I = -2 


L100 - 2(0)] + C[100 - 2(0 )] 3 / 2 => C(100 ) 3 / 2 = -100 

(100 — 2t) 3 / 2 , dy _ r\ dy _ _^ _ (|) (100 - 2t) 1 / 2 (-2) 


= -2+ 3v/l °g~ - = 0 => 20 = 3t/IÓÓ 
=> t « 27.8 min, the timt 
y(27.8) = [100 - 2(27.8)] 


2 t 


10 

400 = 9(100 


dt 

2 t) 


dt 

400 = 900 


18t 


10 

-500= 18t 


t ss 27.8 min, the time to reach the máximum. The máximum amount is then 

[ 100- 2(27.8)] 3 / 2 

10 


14.8 Ib 


16. Let y = y(t) be the amount of carbón monoxide (CO) in the room at time t. The amount of CO entering the 
room is (TOo x tü) ~ Tüüo ft 3 / m in, and the amount of CO leaving the room is ( 4555 ) (y/j) = pTooo ft 3 /min. 


Thus, $ = 12 


dy 


1000 15,000 


+ 


12 


P(t) = 


, Q(t) = 


12 

1000 


15,000 

v(t) = e t/15 ' 000 


15,000 y — 1000 *w - 15,000 
=> y = jzkn fmo et/15 ’ 000 dt =k y = e-‘/«,ooo (12^00 e , / 15 , 00 o + c ) = e -t/i5,ooo (180e ./u,ooo + c) . 

y(0) = 0 0 = 1(180 + C) => C = -180 y = 180 - 180e^ 15 ’ 000 . When the concentration of CO is 0.01% 

in the room, the amount of CO satisfies 45 ) J() = ygg => y = 0.45 ft 3 . When the room contains this amount we 
have 0.45 = 180 - 180e-O 5 ’ 000 => = e -t/i 5 ,ooo ^ t = _i 5> 000 ln (« 37.55 min. 


180 


9.4 GRAPHICAL SOUTIONS OF AUTONOMOUS EQUATIONS 


L y' = (y + 2)(y-3) 

(a) y = —2 is a stable equilibrium valué and y = 3 is an unstable equilibrium. 

(b) y" = (2y - l)y' = 2(y + 2) (y - ±) (y - 3) 


y’> 0 

i y'< 0 

i y'> 0 

i ^ 

é 4 1 1 

X 1 -w . 

-4 

* < 0 ! 2 

T 4 > 5 

y"< 0 

j y*>0 ¡ y'<0 

¡ y"> 0 


0.5 
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(c) 



2 . y' = (y + 2)(y — 2) 

(a) y = — 2 is a stable equilibrium valué and y = 2 is an unstable equilibrium. 

(b) y" = 2yy' = 2(y + 2)y(y — 2) 


y' > 0 

y' < 0 

y' > 0 

-3 * - 

[2 -1 1 

i° i i 

3 * 

O 

V 

V» 

V 

o 

o 

V 

v> 

y" > 0 


(c) 



3. y' = y 3 -y=(y+l)y(y-l) 

(a) y = — 1 and y = 1 is an unstable equilibrium and y = 0 is a stable equilibrium valué. 

(b) y" = (3y 2 - l)y' = 3(y + 1) (y + ^)y(y - (y - 1) 


y' < 0 

y' > 0 

y' < 0 

y' > 0 

-1.5 * J 

y" < 0 

1 -0.5 ' 

y" > 0 

y" < ( 

y" > C 

O 

V 

lO £*■» 
r O . 

1.5 

y" > 0 


1 1 

a /3 y/Z 


y 
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4- y' = y(y-2) 

(a) y = 0 is a stable equilibrium valué and y = 2 is an unstable equilibrium. 

(b) y" = (2y — 2)y' = 2y(y — l)(y — 2) 



(c) 



y > 0» y*> o 
>•' < 0. y* < 0 


y' < 0, 3/ > 0 


/ > 0, j/cO 


5- y'= v/y,y >o 


(a) There are no equilibrium valúes. 



0 12 3 4 

y"> 0 



-2 2 4 6 8 


6 . y , = y- v / y,y>o 

(a) y = 1 is an unstable equilibrium. 

<b> y" = (‘ - 57?) y' = (i - 57?)(y- vy) = (,/y- IKv'y- 0 

y' < 0 j y' < 0 ¡ y' > 0 



0.25 
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(c) 


y 



7. y' = (y- l)(y-2)(y-3) 

(a) y = 1 and y = 3 is an unstable equilibrium and y = 2 is a stable equilibrium valué. 

(b) y" = (3y 2 - 12y + 11)(y - l)(y - 2)(y - 3) = 3(y - l)(y - ^)(y - 2)(y - 


1 

y<o 

• y>o i 

y'< 0 

a 

y>o 

i 

0 


i? ; 2 ? 

! 3? 


4 


y”< 0 

1 y">0 i y*<0 1 

y"> 0 i y"<0 1 

y"> 0 





6 + ^~2.58 




3 3 


(c) 




4 

3.5 



—■— y'< 0, y"< 0 


-----y'<0,y">0 


—■— U y’> 0, > ,r < 0 

——‘■*^1 

—-- y'>0,y'>0 

0.5 

N,^ y'<0,y ff <0 

-1 

1 2 3 ^ 

3 / 

= y - y 

(N 

II 


(a) y = 0 and y = 1 is an unstable equilibrium. 

(b) y" = (3y 2 - 2y)(y 3 - y 2 ) = y 3 (3y - 2)(y - 1) 


• VJ_ 

A 

o 

y' < 0 

y' > 0 

o 

V 

"V. 

O 

A 

V» 

y" < 0 

t 

O 

A 

V> 


2/3 


(c) 


y 
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